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Abstract 

We consider nonresonant contributions in the Dalitz-plot analysis of — > /jvr — > 7r~^TT~iT^ decay and 
their potential impact on the extraction of the Cabibbo-Kobayashi-Maskawa parameter a. In particular, 
we examine the role of the heavy mesons B* and Bq, via the process B — > tt{B*, Bq) — > 7r"'"7r~7r'^, and their 
interference with resonant contributions in the /o-mass region. We discuss the inherent uncertainties and 
suggest that the effects may be substantially smaller than previously indicated. 
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I. INTRODUCTION 



The recent observation of CP violation in the 5-meson system, reahzed through the measurement 
of a nonzero, time-dependent, CP-violating asymmetry in the process P*^(P°) — > J/ipKg (and 
related ones) [1], heralds a new era of discovery. The result yields a value of sin(2/3) in accord 
with standard model (SM) expectations [2], where P, defined by exp(i/9) = —V*^V^^/ (V^IV^^) , is an 
angle of the unitarity triangle, V^j being an element of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix [3]. Ascertaining the presence of physics beyond the SM thus demands the determination 
of all the angles of the unitarity triangle. 

In this paper, we consider the decays P°(P°) ^ pvr — > 7r"'"7r~7r°, as a Dalitz-plot analysis of the 
possible pn final states, under the assumption of isospin symmetry, permits the determination of 
the CKM parameter a [4], where a — n — P — ^ and exp(i7) = —V*jy^^/{V*jy^^). Our interest 
is in assessing the size of the nonresonant contributions which could possibly obscure the analysis, 
and in ameliorating their impact. Indeed, the strategy for the extraction of a relies, in part, on 
the assumption that the p mesons dominate the Svr final state. There are, however, empirical 
indications that this assumption may not always be warranted. For example, combining the CLEO 
measurements of the branching fractions, B{B'^ p'^n^) = {27.6^.^1 ±4.2) x 10~^ and B{B~ — > 
p°7r-) = {WAtli ± 2.1) X 10-6 [5], with the BABAR result B(P° ^ pV) = (28.9 ± 5.4 ± 4.3) x 
10-6 [6] yields 

where we have added the errors in quadrature and ignored correlations. These ratios are small [7] 
with respect to simple theoretical estimates, which give 7^ ~ 6 [8]. An interesting possibility for 
the resolution of this discrepancy has been suggested in Refs. [9, 10], whose authors investigate 
the possible backgrounds to P — > p7r — > Stt decay which arise from contributions mediated by 
other resonances. They find that the light a resonance, a broad / = J = enhancement in tttt 
scattering, as well as the heavy-meson resonances B* {J^ = 1~) and Pq {J^ = O"*"), can modify the 
P — s> Stt branching ratios in the p-mass region and give rise to values of TZ crudely compatible with 
the empirical value of Eq. (1), given its large error. In particular, the contribution of B~ — > a7r~ 
decay significantly enhances the effective B~ — > p°7r- branching ratio and lowers the value of TZ. 
Analogously, the a modestly impacts the P° — > p°7r° branching ratio [11]; let us consider the issues. 

The analysis of P°(P°) — > p7r — > n'^n'Tr^ decay posits a two-step process, that is, that the 
amplitude for P° 7r''"7r-7r° decay can be written as 

A{B' ^ 7r+7r-7r°) = /+a+_ + /_a_+ + f,a,, , (2) 

where a- = A{B^ — > pV-') and is the vector form- factor describing tttt [4]. An analogous 

construct can be made for B^{B^) air — > tt+tt^tt^ decay, which contains the scalar form-factor 
describing a — > ti^ti' . It is evident that the manner in which the a populates the p phase-space will 
depend on the amplitude for P° — > (77r° decay, as well as on the accompanying scalar form-factor. 
The (7 is a state of definite CP, so that the isospin analysis of Ref. [4] can be enlarged to include 
it [11]; nevertheless, the analysis relies on the form factors adopted for the p* — > tttt and cr — > tttt 
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processes. The resonances of interest are broad, so that Breit-Wigner form-factors arc gcncraUy 
insufficient: they do not satisfy general theoretical constraints, such as analyticity and unitarity, 
over the tttt invariant-mass interval needed. As discussed in detail in Ref. [11], the differences are 
striking for the scalar form-factor, and the resulting numerical impact on S — > Stt decay is sizable. 
In contrast, the numerical differences for the vector form-factor are not large. 

The purpose of this paper is to extend the work of Ref. [11], which deals exclusively with the 
p and a contributions. We incorporate the B* and Bq contributions suggested in Ref. [9], as the 
effects they find in the B^ p^n^ channel are considerable. In this paper, however, we show that 
the off-shell nature of the B* and Bq weak and strong vertices adds considerably to the uncertainty 
of the estimate of Ref. [9] and may well reduce these contributions significantly. Nevertheless, we 
also explore kinematical cuts which would be useful in reducing the impact of these effects in the 
p-mass region. 

We begin in Sec. II with the weak, effective Hamiltonian and the matrix elements pertinent to 
our calculations. Subsequently, in Sec. Ill, we derive the amplitudes associated with the various 
contributions of interest in the p-mass region oi B ^ decay. We discuss our numerical results 
in Sec. IV and conclude in Sec. V. 



II. EFFECTIVE HAMILTONIAN AND MATRIX ELEMENTS 



The effective, \AB\ — 1 Hamiltonian for b dqq decay is given by [12] 

10 



^eff 



G 



F 

a/2 



(3) 



where Gp is the Fermi coupling constant, = V^^y*^ are CKM factors, Q are Wilson coefficients, 
and Oj are four-quark operators. The expressions for Q and are detailed in Ref. [12], though 
we interchange C]0\ <-> C2O2, so that ~ 1 and > C2 . We neglect the electroweak- 
penguin operators 0-j ... because their coefficients C7 ... arc smaller than the others. In the decay 
amplitudes that we derive, the enter through the combinations = C^-\- C^_^i/N^ if i is odd 
and = + C^_]^/N^ if i is even, where iV^ = 3 is the number of colors. 

The diagrams contributing to the B ^ "in amplitudes considered here, as shown in Fig. 1, each 
have a strong vertex and a weak vertex, where the latter describes the transition — >■ M;^M2, 
in which is a heavy meson containing a h quark and 2 s-^® light mesons. The amplitude 
corresponding to the weak vertex is given by 

A{M, ^ M,M,) = {M,M,\H,^\M,) . (4) 

To evaluate this, we adopt the naive factorization approximation, following earlier calculations [9, 
10, 11] to which we compare. 

The relevant matrix elements are 

(7r-(p)|dVLii|0) = V2{'k\p)\u^^Lu\Q) = i/^p'^ , 

^ (5) 
{p-{p,e)\dYu\Q) = V2{p\p,e)\u^^u\Q) = , 
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B a- < B 

TT 

(a) (b) 
FIG. 1: Diagrams contributing to ^ Svr, decay with each square denoting a weak vertex. 

g^(p+(p,e)|n7^L6|5°(A:)) = -2i<^^(g2) ■ g , 

Ml -Ml 



{p)\urLb\B'{k)) = (k+pfFf-^iq') + ^ ^ q''{Fr^{q^)-F^--{q')) , (6) 

qA 

q^{a{p)\d^^Lh\B\k)) = -i (M| - M^) Fo^-(g2) , 



q'^{n+{p)\u^^Lb\B*'{k,ej,,)) = 2iV2 A^' ^^q') M e ■ q , 
q^{'K^{p)\uYLh\Bl{k)) = -i (m|^ - M^) F^^^-{q') , 

where and are the usual decay constants, q = k — p, and L = 1 — 75. The various A^i^q"^) and 
-^01(5^) form factors. Other meson-to-meson matrix elements can be determined using isospin 
symmetry. In our phase convention, the meson flavor wave functions are given by vr^ = ud, \pi -k^ = 
uu — dd, 7T~ = du, B^ = bd, B~ = bu, and similarly for the p, B*, and B^. This implies that we 
have, for example, {7i+\u-f^b\B^) = -^2 {n'^\dj^b\B^) = +V2 {7r^\uj^b\B-) = {7i-\d-f^b\B'). We 
now employ these matrix elements to realize amplitudes for B ^ Sir decays. 



III. AMPLITUDES 

Practical considerations drive our interest in the 7T~^7[~7r^ and tt^tt^vt^ decay modes; we shall 
not consider the Tr'^Tr^vr^ ones. We write the amplitude for B^ 7r+(p^) TT^{p_) tt^{Pq) decay as a 
coherent sum of the p, a, B*, and Bq amplitudes, namely, 

A^-' = a;-' + At-' + A+r° + . (8) 

For B^ n^^Pi) 7r^(p2) '^^{p+), the amplitude A ^ can be constructed in an analogous manner. 

We consider first the B ^ pn ^ contributions, represented by the diagram denoted by 
"p" in Fig. 1(a). For each p* diagram and 3tt state, the amplitude is written as a product of an 
amplitude for the B pV-' weak transition and a vertex function F^^^ describing the p* tttt 
form factor. Were the p a narrow resonance, the Breit-Wigner (BW) form 



would suffice, where ^/s is the invariant mass of the 27r system and Qp is the p ^ tttt coupling 
constant. However, since the p is not narrow — its width is some 20% of its mass — this form 
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must be generalized to accommodate known theoretical constraints over the region in s for which 
it is appreciable. For example, unitarity and time-reversal invariance compel the phase of Tp^^^s) 
to be that of L — 1, I — 1 tttt scattering for s < (M^ + M^ f, for which the scattering is elastic. 
Moreover, the imaginary part of r^^(s) must vanish below physical threshold, s — 4M^. For a 
detailed discussion with references to earher work, see Refs. [11, 13]. Following Ref. [13], we have 

-F (s) 

rp..(«) = , (10) 

J fry 

where Fp{s) is the vector form-factor of the pion and /^^ is the p-7 coupling constant. The pa- 
rameters in -Fp(s) are determined by fitting to e+e^ tt+tt^ data; what is important is that the 
parametrization itself is consistent with theoretical constraints. The value of /^^ is determined 
from the p — > e'^e~ width, which, in turn, is extracted from the e'^e" — > tt'^tt" cross section 
at s = Mp [13, 14]. The overall sign is chosen so that Eq. (10) is equivalent to the BW form, 
Eq. (9), as s — > . At s = the BW form is compatible with the various theoretical con- 
straints. In our numerical analysis, we adopt the "solution B" fit of Ref. [13] for F^, for which 
fp^ = 0.122 ± 0.001 GeV^ [14]. Alternatively, a BW form with a running width rp(s), chosen to 
be compatible with the form of the tttt phase shift (in the crossed channel) as s — > 4M^, is given 
in Ref. [15]. However, the numerical differences between this form and the one we have chosen are 
small [11]. 

For the decay amplitudes, after summing over the p polarizations, we find 

"° = V~ {s+- - s+o) ^pnnis-o) + (s-o - «+-) ^p^wis+o) - {s_Q - s+o) Tp^^{s+_) , 

r 1 ^^^^ 

where s^i = {Pk+pf, with 

V- = ^(A„ai-A,aJ/,Ff-, r)+ = ^ [A„a, - A, (a^ - a^i?,)] /.M^^f , 

= ^{[X^a, + \{a,-a,R^)]f,MpAl^'^ + {X^a, + \a,) fpFf^] , (12) 

[K «i - \ («4 - «6 R,)] LMp Al' + (A„ + A, a J F^^ } . 

Here = A^^^M^) and Ff^ = Ff^^(M2), whereas = M^/[{mf, + m)m] — note that 
we work in the isospin-symmetric limit, for which rfi = m„ = . The relative signs between the 
different terms in Eq. (11) follow from the pnn couplings^ 

(7r°(po)7r^(p±)|p^) = ^QpSp ■ {p±-Po) , 

(13) 

(7r+(p+)7r (p_)|p°) = gpep-{p_-p^) , 



^ We use the notation {M2M^\M-i) = {M^M^\Us^^o^^\M-^) . 
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which follow, in turn, from the phase conventions we have chosen for the flavor wave functions: 
Itt^) = =f|-^ = 1,-^3 = ±1) and |7r°) = \I = 1, = 0), and similarly for the p states. Our 
amphtudes agree with those of earher calculations [9, 11, 16]. 

We turn next to the a "meson" contributions, represented by the diagram denoted by "cr" in 
Fig. 1(a). We use the a to denote a two-pion state with total isospin 7 = and total angular- 
momentum J = ; it need not be a "pre-existing" resonance, but, rather, can be generated dy- 
namically by the strong pionic final-state interactions in this channel [17]. The peak of the broad 
enhancement associated with the a is close to the p in mass, so that the decay B air Stt can 
populate the B ^ pir phase space [10]. As in the p case, the amplitudes for 5 — > cttt — > Stt decays 
are written as a product of an amplitude for the 5 — > cttt weak transition and a vertex function 
r^^^ describing the cr — > tttt form factor. We write [11] 



where is defined as 



r...(5) = xrris) , (14) 



(0|Jd|7r+(p+)7r-(p_)) = ^JhUs+-)B, . (15) 



We note that B^ = M^/(2m) is the vacuum quark condensate and x is a normalization constant, 
to be discussed shortly. For our numerical work in the next section, we adopt the r"(s) as derived 
in Ref. [18], after Refs. [17, 19, 20]. The calculated form factor is reahzed in a chiral, unitarized, 
coupled-channel approach; at low energies, the form factor is matched to the one-loop-order expres- 
sion in chiral perturbation theory [18, 21]. The resulting form factor is consistent with low-energy 
constraints and is comparable to the scalar form-factor which emerges from the dispersion analysis 
of Ref. [22]; however, it is notably different from the Breit-Wigner form adopted in Refs. [10, 23] to 
study the role of the a in B and D decays into the Stt final state. That is. 



X 



r^^(s) = r (s) = M^J^^^ML (16) 

where the coupling g^^^ = (7r+7r~|cr) is determined from the o" — > tttt decay rate. For S — > Stt 
decay, the numerical changes arising from the use of r^^^(s) in place of the BW expression are 
significant [11], as we will see here as well. We determine the normalization x by requiring that [11] 

which equates [Fo-^ttCs)! to its BW counterpart at s = M^. The values of Mg. and F^. are extracted 
from fits of F^^(s) to D — > Stt decays [23]. The normalization condition is motivated by noting 
that the modulus of F"(s) is peaked near s — M^, whereas the normalization of T^{s) is sensitive 
to the values of certain, poorly known low-energy constants [11]. We emphasize that and F^^ 
appear merely in the normalization of F^^^. 

The resulting decay amplitudes are then 

Ar' = r_(«+-) , = (r...(«i+) + r_(«2-K)) , (is) 



6 



where 

[K «2 + \ («4 - «6 R,)] {Ml - Ml) f^Ft - A, a, (Ml - Ml) 

(19) 

t-^[ [K a, - A, (a, - a, R^)] (M| - M,^) f^F^ + A, (M| - M,^) Fo^^} , 

with F(f'^ = Fo^^^(M2) and F^" = F^^^^Ml). From Eqs. (14) and (15), it follows that 
{a\dd\Q) = Ml / x'lTT') ■ We agree with the weak amplitudes of Ref. [11], but disagree with those 
of Ref. [10] in that our r]^ and fj^, neglecting penguin terms, are smaller and larger, respectively, 
than theirs by a factor of ^/2. 

We now evaluate the B* and Bq contributions, whose diagrams are shown in Fig. 1(b); we 
suppose that other excited i?-meson states could also contribute, but we expect that their larger 
masses ought to make them less important [9]. Presently, no reliable data exist on the widths of 
these heavy mesons, so that their values have to be calculated. Recent estimates [24, 25] suggest that 
the B* is a very narrow resonance, whereas the Bq is less so, its width being some 6% of its mass. 
Nevertheless, the resonances are sufficiently narrow that it is reasonable to adopt a Breit-Wigner 
representation for the propagators of these mesons, as in Ref. [9]. In the combined heavy-quark and 
chiral limit [26], the strong couplings connecting the (B*,Bq), B, and tt mesons are [9, 24, 27]^ 

{B-{p')n^p)\B*%k,s)) = J_9^/M^e.p , (20) 

/ IT 



h M^M^ f.2_M2 

Using isospin symmetry, we derive 

(S-7r+|5*°) = -\/2(^V|B*°) = {B\-\B*-) (22) 
and analogous relations for (SttISq). We then obtain 

n(s_o, s^_) + n(s_o, s+o) n(s+_, «_„) 



4+-0 



(23) 



At 



-0 



si+ - Ml, + iP^.M^, - Ml, + iP^. Mj,, ' 

+ ir- \ {mi, - Ml 



' ^-0 - Ml^ + iP^^M^^ - Ml^ + iP^^M^^ j V2 

A'" A'" 



(24) 



I — + 



(m^ - Ml) , 



We note that {B*- {k,e)-K+\B^) = -{B--K+\B*^{k,e)) and {Bo{k)-K+\B^) = -{B--K+\B^{k)). 
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where 



K = -4Gp [\a, - \ (a^ - a^R^)] gM,,.^ M^M^, A^'^ , 



V Bg 



(25) 



(26) 



and the sum over B* polarizations yields 

, , (Ml- Ml-u)u ^ V 

Note that ^ = -^^(M^) and F(f"" = F(f°^"(M2). Our expressions for in Eqs. (23) 

and (24) disagree with those in Ref. [9] in that the factors of l/-\/2 are missing in their formulas, 
and that the minus sign in the middle of the big brackets in Eq. (24) is opposite to theirs. However, 
our expressions for A in Eqs. (23) and (24) agree with theirs. 



IV. NUMERICAL RESULTS AND DISCUSSION 

We begin by listing the parameters that we use; we conform with the parameter choices of 
Rcfs. [9, 11], in order to realize a crisp comparison with their results. In specific, the Wilson 
coefficients we use are 

Ci = 1.100, C2 = -0.226, C3 = 0.012, C4 = -0.029, C5 = 0.009, Cg = -0.033. (28) 

For the CKM factors, we adopt the Wolfenstein parametrization [28], retaining terms of O(A^) in 
the real part and of O(A^) in the imaginary part, to wit, 

K<i = I-AV2 , = ^A^[p-i77(l-AV2)] , T4d = AA^(l-p-i77) , V,, ^ I , (29) 
and using 

A = 0.2196 , p = 0.05 , 77 = 0.36 , A = 0.806 . (30) 

For decay constants, light meson masses, and resonance parameters, we have 

f^/V2 = 92.4 MeV , = 139.57 MeV , 

fp = 0.15 GeV^ , Mp = 769.3 MeV , = 150 MeV , Qp = 5.8 , (31) 

= 478 MeV , F^ = 324 MeV , g^^^ = 2.52 GeV . 
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The decay constants and are associated with tt^ and decay, respectively. We neglect 
isospin-violating effects throughout, so that M^± = M^o = M^, M^± = M^o = M^, as well as 
Afgo = = M^. Moreover, m — 6MeV. The B* and B are degenerate in the heavy-quark 

hmit, so that we neglect their mass difference as well. We also neglect the hfetime difference between 
the S° and B~, setting r^o = Tg- = t^- For the B and related mesons, we have 



= 5.279 GeV, = 1.6 x IQ-'^s , m,, = 4.6GeV , 

r^. = 0.2 keV , = 5.697 GeV , T^^ = 0.36 GeV , 



(32) 



and use 



^ = 0.6 , h ^ -0.7 . (33) 

The heavy-to-light transition form factors are given by 



= 0.29 , F^"" = 0.37 , Ff"" = 0.37 , F^" = 0.46 , 
= 0.16 , Fo'"" = -0.19 . 



(34) 



Finally, for the vector and scalar form-factors, r^^^(s) and r^^^(s), respectively, we follow the 
treatment of Ref. [11]. The -Fp(s) parametrization we adopt was fit to e+e~ tt+tt" data in the 
elastic region [13], 2M^ < ^/s < M^ + M^, only, so that for larger values of s wc use a Breit-Wigner 
form, matched to the value of Tp^^{s) at ^/s = + M^. That is, for ^/s > 923 MeV we employ 
rp..(«) = k(M2 -s)+ iqr^M^]5^/[(M2 - sf + TlMli with 0.929 and c, ^ 1.29. For 

the scalar form- factor, wc employ the r"(s) derived in Ref. [18], which is valid for ^Js < 1.2 GeV. 
The normalization of Eq. (17) implies that x = 20.0 GeV~^. For ^/s > 1.2 GeV, we match to the 
asymptotic form of T^^^[s) [22], as detailed in Ref. [11]. 

To obtain branching ratios for i? ^ 37r decay in the p-mass region, we integrate over the region 
of phase space satisfying the requirement that two of the three pions reconstruct the p mass within 
an interval of 25, as was done in Refs. [9, 11]. This amounts in each case to calculating the effective 
width 

reff(^^p(Pi+P2)7r(p3)) = r(5^7r(pi)7r(p2)7r(P3)) (35) 

We choose 5 = 0.3 GeV, following earlier work [9, 11]. 

For crisp comparison with Ref. [9], we begin by computing the effective branching ratios arising 
from the use of Breit-Wigner forms, as in Eqs. (9) and (16) for the p and a, respectively, throughout. 
The various contributions, reflective of the enumerated terms in Eq. (8), are reported in Table I. 
There are differences between our results for the p, p+B* , and p+B*+Bq contributions and the 
corresponding ones in Ref. [9]. The differences are, however, not large and arise in part from 
missing factors in the formulas for the B* and Bq amplitudes, which we delineated in the last 
section. In contrast, as pointed out in Ref. [11], the a effect on the B decay is much bigger than 
that found in Ref. [10], because our a amplitude is larger than theirs by a factor of This is 
evident in the p+a and p+a+B* columns. Our results agree with those Ref. [11], to the extent 
that they are apphcable; we note that Ref. [11] neglects penguin contributions altogether and deals 
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TABLE I: Effective branching ratios for B ^ pix decays, as per Eq. (35), with 5 = 0.3 GeV. Breit-Wigner 
form factors are used throughout, noting Eqs. (9) and (16) for the p and a contributions, respectively. Ah 
branching ratios are reported in units of 10~^. 



Decay mode 


P 


a 


B* 


Bo 


p + B* 


p+B*+B„ 


p + a 


p + a + B* 


p + (7 + B*+Bo 




16.0 


0.0003 


0.54 


0.009 


16.5 


16.3 


16.0 


16.4 


16.3 




4.76 


0.0003 


0.13 


0.020 


4.98 


4.98 


4.78 


5.00 


5.00 




0.91 


0.045 


0.39 


0.016 


1.43 


1.29 


0.93 


1.59 


1.43 


B-^p°Tr- 


4.10 


5.18 


2.71 


0.107 


7.42 


8.45 


8.83 


7.67 


7.92 


■R 


5.1 








2.9 


2.5 


2.3 


2.8 


2.7 



exclusively with the p and a contributions. The last column of Table I contains the sum of all the 
contributions, p+a+B*+BQ . Overall, it is apparent that the effect of the Bq is smaller than that 
of the other contributions, although it is not negligible. Finally, in the last row, we collect the ratios 
of branching ratios TZ defined in Eq. (1). These results show that the inclusion of the a and B*, 
either individually or together, makes the estimated value of TZ consistent with the empirical one, 
given its large error. 

We now proceed to compute the effective branching ratios with the p and a form-factors, Eqs. (10) 
and (14), which we advocate. These results are presented in Table II. The results without the a 
contributions change httle, as the vector form-factor is not terribly different from its BW counter- 
part [11]. In the presence of the a, this similarity persists for the S° decays, but, in contrast, the 
B~ branching ratios are significantly increased compared to the corresponding ones in Table I. This 
effect also tends to diminish the relative impact of the B* and Bq contributions on the p^n^ mode, 
though the heavy mesons persist in making a substantial impact on the effective branching ratio 
for the p^n^ mode. 



TABLE II: Effective branching ratios for B ^ pn decays, as per Eq. (35), with 6 = 0.3 GeV. We adopt 
the p and a form factors, Eqs. (10) and (14), respectively, which we have advocated. All branching ratios 
are reported in units of 10~^. 



Decay mode 


P 


a 


B* 


Bo 


P + B* 


P + B*+Bq 


P+CF 


p + a + B* 


p + a + B* +Bo 


B°^p-7:+ 


16.0 


0.001 


0.54 


0.009 


16.6 


16.4 


15.9 


16.5 


16.3 




4.76 


0.001 


0.13 


0.020 


4.90 


4.93 


4.80 


4.94 


4.98 




0.86 


0.065 


0.39 


0.016 


1.35 


1.21 


0.91 


1.47 


1.33 


B-^p^TT' 


4.06 


7.66 


2.71 


0.107 


7.20 


8.25 


11.1 


11.9 


12.7 


n 


5.1 








3.0 


2.6 


1.9 


1.8 


1.7 
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Were the heavy-meson contributions to the p^n^ mode seen in Table II as large as we have 
estimated, the impact on the Dalitz-plot analysis to extract a from B^{B^) vr+vr^Tr^ decays would 
be significant [9]. Since the B* and Bq masses lie outside the phase-space region of 5 — * Svr, their 
effects behave as part of the nonresonant background, but are not uniform and obviously interfere 
with other contributions. The manner in which the contributions are distributed throughout the 
Dalitz plot is shown for B^ — > 7r+7r~7r° decay in Fig. 2; the heavy-meson contributions preferentially 
populate the edges of the Dalitz plot, in which the p contributions lie as well. In B~ —>■ tt'^tt'tt' 
decay, the distribution of the heavy-meson contributions is somewhat more uniform, as illustrated 
in Fig. 3. 




FIG. 2: The B* and Bo contributions to 5° 7r+(p+)7r- (p_)7r°(po) decay, specifically, \A^7'^ + A 
(in dimensionless units) as a function of its arguments s_,_q and s_q, both in units of GeV^. 
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FIG. 3: The i?* and i?o contributions to i? ^ vr (pi)vr {p2)7T {p+) decay, specifically, 
(in dimensionless units) as a function of its arguments si+ and S2+, both in units of GeV^. 
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We now proceed to consider the reliability of the estimates we have effected. Let us first note 
that the parameters g and h of the strong heavy-meson couplings in Eqs. (20), (21) assume the 
values given in Eq. (33) — these reflect the upper limits of their estimated ranges [9, 24].^ Thus, the 
results we find with these parameters can be regarded as extremal estimates (although variations 
in other numerical inputs, such as the form factors, do introduce further uncertainties). Choosing 
central values of g and h in their estimated ranges decreases the heavy-meson effects by up to some 
50% [9], as explicitly shown in Table III. 

TABLE III: Effective branching ratios for B ^ pir decays, as in Table II, except that g = 0.40 and 
h = —0.54 have been used. 



Decay mode 


p 


B* 


So 


p + B* 


p + B*+Bo 


P + (T 


p + a + B* 


p + a + B*+Bo 




16.0 


0.24 


0.005 


16.3 


16.1 


15.9 


16.2 


16.1 




4.76 


0.06 


0.012 


4.82 


4.87 


4.80 


4.86 


4.91 




0.86 


0.17 


0.009 


1.10 


1.03 


0.91 


1.20 


1.13 




4.0(i 


1.20 


().0()4 


5..-).-) 


(i.l2 


11.1 


11.0 


11.4 


n 


5.1 






3.8 


3.4 


1.9 


1.9 


1.8 



Moreover, the relative signs chosen for the p, cr, and heavy-meson contributions will impact the 
numerical values of the effective branching ratios. As noted by Ref. [9], the relative sign of the 
B* and Bq contributions is fixed in the heavy-quark and chiral hmits. The relative signs of the 
heavy-meson, p, and a contributions, however, are less clear. We define the p ^ tttt coupling as 
per Eq. (13), after Ref. [9, 11], though we note that a chiral Lagrangian analysis suggests that 
the relations of Eq. (13) should possess an additional overall sign. With this modification, the 
branching ratios for the p+B*+Bq combination in Table II typically become smaller by no more 
than 15%. However, the p+a results in B^ — >• p°7r° and B~ — >• p°7r~ become some 3% and 10% 
larger, respectively. The impact on the p+a+B*+BQ results is mixed, leading to a suppression of 
about 10% in the p°7r° mode and an enhancement of 2% in the p°7r~ mode. 

Kinematical cuts can mitigate the impact of the heavy-meson and a contributions. Since the 
p^n^ modes are little affected by these notions, we evaluate only the p°7r° and p°7r~ modes. We 
try two different sets of kinematical cuts. For the first one, we set S — 0.15 GeV = L^ and report 
our results in Table IV. The relative suppression of the heavy-meson and a contributions is quite 
modest, if it exists at all. For the second set, we impose not only a d cut but also a cut on cos^, 
where 6 is the helicity angle, defined as the angle between the direction of one member of a pion pair 
from p decay and the direction of the parent S-meson evaluated in the pair's rest-frame. Since the p 
contribution has a cos^^ distribution in B ^ tt+tt^tt^ decay [4], larger values of | cos^j enhance the 
p contribution. Interference effects in the B" — > 7r"^7r~7r~ channel will make this cut less effective. 



We also note, however, that the g vahic in Eq. (33) is, by virtue of heavy-quark symmetry, favored by the recent 
measurement of the D* Dtt width [29], which yields g = 0.59 ± 0.01 ± 0.07. 
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TABLE IV: Effective branching ratios for B ^ pir decays, as in Table II, except that S = 0.15 GeV has 
been used. 



Decay mode 


P 


a 


B* 


p + B* 


p + B*+Bo 


p + a 


p + a + B* 


p + a + B*+Bo 




0.33 


0.029 


0.20 


0.55 


0.49 


0.36 


0.59 


0.53 




3.3G 


3.46 


L38 


4.82 


5.39 


6.46 


7.56 


8.17 



TABLE V: Effective branching ratios for B ^ pn decays, as in Table II, but with the additional kinematical 
cut I cos^l > 0.4 as explained in the text. 



Decay mode 


P 


a 


B* 


p + B* 


p + B* +Bq 


p + o 


p + a + B* 


p + a + B* +Bo 




0.84 


0.039 


0.27 


1.22 


1.11 


0.87 


1.29 


1.18 




3.81 


4.79 


1.71 


5.92 


6.58 


7.97 


8.63 


9.13 



We set 6 = 0.3 GeV and | cos^^| > 0.4, and collect the results in Table V. Comparing to Table II, 
the 6 cut is seen to decrease the relative size of the a background, as discussed in Ref. [11]. The 
helicity-angle cut only modestly reduces the a contribution in — p^ir^ decay; however, an 
assumption of p dominance is only needed when one employs an isospin analysis to extract sin(2Q;) 
from — > 7r''"7r~7r° decay, as detailed in Ref. [4]. For the p°7r° mode, the chosen cut does 
significantly reduce an already small contribution. Were the a contribution to the p^n'^ mode much 
larger than we estimate, then a full partial- wave analysis to separate the s- andp-wave contributions 
could be both practicable and necessary. 

Finally, we must discuss a tacit assumption we have made in the estimation of the B* and Bq 
contributions, which is made in Ref. [9] as well. That is, in realizing the diagrams of Fig. 1(b), we 
have treated the strong {B*, Bq)Bt: and weak [B* , Bq) — > tttt vertices as if the B* and Bq mesons 
were on their mass shell. This assumption is compatible with the assumed use of the combined 
heavy-quark and chiral limits in the treatment of the strong {B* , Bq)B'k vertices. However, neither 
assumption is appropriate for B ^ pir decay. That is, for the ^ 37r decays of interest, we 
require that two of the three pions have an invariant mass comparable to that of the p meson. 
This implies that in most of the relevant phase-space region the mediating heavy-mesons carry s 
values much smaller than their squared masses — they are highly off-mass-shell. Moreover, the 
bachelor tt is never soft in this kinematical region. Thus the combined heavy-quark and chiral 
limits are used beyond their range of validity. These effects modify the vertices we have assumed 
in Eqs. (20), (21) and Eq. (7). Unfortunately, the needed off-shell extrapolations cannot be done 
reliably, although we would generically expect this effect to suppress the numerical importance of 
the B* and Bq contributions. For example, the form factors of Eqs. (7) now depend on both (f' 
and /c^; the vertices are only "half" off-shell, so that does not enter, as the final-state tt is on 
its mass shell. Moreover, additional form factors appear. To illustrate, we note that the general 
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parametrization 



1^ 



(36) 



predicated by an assumption of Lorentz invariance yields 

ct(,^'^{v)\u^^Lh\Bl{k)) = -i (m|^ - Ut) F^^^\q\ k') + i (m|^ - k') Pf^^'^iq', k') (37) 

for the half-off-shell matrix element of interest. The matrix element is a linear combination of 
signed, uncertain contributions, so that its sign is ultimately unclear. Similar considerations apply 
to the S* — > TT matrix element, as well as to the strong vertices of Eqs. (20), (21). In the treatment 
of Ref. [30], an off-shell extrapolation of Eq. (20), in the kinematic region of interest, is effected 
through the replacement \/M^M^ — > -^/MB^fs- To assess the impact of these considerations on 
the numerical results we have reported, we shall adopt a similarly ad hoc prescription. Thus, we 
perform the replacement 



m; 



3/2 



.3/4 



(38) 



in the numerator of the B* amplitudes in Eq. (23), so that the "off-shellness" of both the strong and 
weak vertices is taken into account. We neglect the Eq in this simple numerical estimate, as its effect 
was rather small to start with. We calculate the corresponding branching ratios and collect the 
results in Table VI. Our simple prescription leads to a dramatic reduction of the B* contributions, 
as a comparison with Table II makes clear. Note that the computed values of TZ are still consistent 
with the empirical ones, as a reduction in TZ is still realized through the a contributions. Although 
we cannot draw firm conclusions from this simple exercise, it serves to illustrate that neglecting the 
off-shell nature of the heavy-meson vertices in the kinematic region of interest could easily lead to 
a considerable overestimate of their effects. 



TABLE VI: Effective branching ratios for B ^ pir decays, as in Table II, except that the off-shellness of 
the B* meson is included as explained in the text. 



Decay mode 


P 


a 


B* 


p + B* 


P+ (7 


p + a + B* 




16.0 


0.001 


0.03 


16.0 


15.9 


16.0 




4.76 


0.001 


0.01 


4.85 


4.80 


4.88 




0.86 


0.065 


0.02 


0.88 


0.91 


0.95 


B-^p'^TT- 


4.06 


7.66 


0.25 


4.43 


11.1 


10.7 




5.1 






4.7 


1.9 


1.9 
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V. CONCLUSIONS 



We have examined resonant and nonresonant backgrounds to i? — > pyr — > Stt decays which can 
potentially impact the extraction of a from a Dalitz-plot analysis of S — > 7r'^7r~7r° decays [4] , as well 
as the value of the ratio of branching ratios we term TZ, as defined in Eq. (1). In particular, we have 
evaluated the effects of nonresonant contributions mediated by the heavy mesons B* and Sq, as 
well as the contributions from the light a resonance via B air Sn decay, in the p-mass region. 
In this, our analysis parallels that of Refs. [9, 10], though it differs fundamentally in two points. 
Firstly, we use the vector and scalar form-factors of Ref. [11], which are consistent with low-energy 
theoretical constraints and thus are suitable for the description of broad resonant structures such as 
the p and the a. The scalar form factor, in particular, is quite different from the Breit-Wigner form 
adopted in other analyses [10, 23] and leads to differing results [11]. Secondly, in the kinematics 
of interest, the B* and are highly off-mass-shell, impacting the strong and weak vertices which 
mediate the B {B*, Bq)t{ —>■ tttttt decay. We find that these effects can reduce the heavy-meson 
contributions substantially. 

Our numerical results show, were we to neglect the off-shell effects we have mentioned, that 
the B p'^TT^ decay modes are little affected by the a and heavy-meson backgrounds, whereas 
the p^n^ mode receives large contributions from the latter. In contrast, the B~ 7r+7r~7r~ decay 
mode contains large contributions from both the a and B*, though the a contributions numerically 
dominate. Effecting a simple model of off-shell effects, we find that the B* effects are substantially 
reduced. The off-shell extrapolation of interest cannot be effected with certainty; nevertheless, our 
estimates indicate that the neglect of this effect may lead to a substantial overestimate of the B* 
contributions in 5 — > Stt decay. The role of the a in lowering the theoretical value of TZ and yielding 
a favorable comparison with experiment persists despite these considerations. 

Note added. Since the submission of this paper for publication, a report by the BABAR Col- 
laboration has appeared [31], giving the experimental bound B{B^ tt+tt^tt^) < 15 x 10^^ at 
90% C.L. This can be used to constrain the contribution of the B*- and fig-pole diagrams. Using 
the g and h values as in Table II, we find B{B~ — > 7r"'"7r~7r~) = 24.8 x 10~^ for the combined 
p+a+B*+BQ contribution, where we have integrated over all the allowed phase-space. Were we to 
use the intermediate values of g and h given in Table III, though such a is not favored by data [29], 
we would obtain B{B~ 7r~^7r~n~) — 18.7 x 10~^. If we use our off-shell extrapolation (neglecting 
the small Bq contribution) and the parameters of Table II, we find B{B~ ti^ti~ti~) = 15.4 x 10^^. 
This comparison supports our assertion: the treatment of the B* vertices in Ref. [9] tends to yield 
an overestimate of their contribution to fi — > Stt decay. On a related note, the failure to confront 
the empirical bound on B{B~ — > K'^K~7r~) decay has been described in recent work by Cheng 
and Yang [32]. 

Acknowledgments We thank Ulf-G. Meifiner for helpful discussions and J. A. Oiler for the use of 
his scalar form factor program. The work of J.T. is supported by the U.S. Department of Energy 
under contract DE-FG01-00ER45832. S.G. thanks the SLAC Theory Group for hospitality and is 
supported by the U.S. Department of Energy under contracts DE-FG02-96ER40989 and DE-AC03- 
76SF00515. 



15 



[1] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 87, 091801 (2001); K. Abe et al. [Belle 
Collaboration], ibid. 87, 091802 (2001). 

See, e.g., Y. Nir, arXiv:hep-ph/0109090; A.J. Buras, arXiv:hep-ph/0109197. 

N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531; M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 
(1973) 652. 

H.J. Lipkin, Y. Nir, H.R. Quinn, and A. Snyder, Phys. Rev. D 44, 1454 (1991); A.E. Snyder and 
H.R. Quinn, ibid. 48, 2139 (1993); H.R. Quinn and J.P. Silva, ibid. 62, 054002 (2000). 
CP. Jessop et al. [CLEO Collaboration], Phys. Rev. Lett. 85, 2881 (2000). 
B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0107058. 

Y. Gao and F. Wurthwein [CLEO Collaboration], arXiv:hep-ex/9904008. 
M. Bauer, B. Stech, and M. Wirbel, Z. Phys. C 34, 103 (1987). 

A. Deandrea et al, Phys. Rev. D 62, 036001 (2000); A. Deandrea, arXiv:hep-ph/0005014. 
A. Deandrea and A.D. Polosa, Phys. Rev. Lett. 86, 216 (2001). 
S. Gardner and U.-G. Meifiner, hep-ph/0112281, to appear in Phys. Rev. D. 
G. Buchalla, A.J. Buras, and M.E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
S. Gardner and H.B. O'ConneU, Phys. Rev. D 59, 076002 (1999). 

S. Gardner and H.B. O'ConneU, Phys. Rev. D 57, 2716 (1998); ibid. 62, 019903(E) (1999). 
P.F. Harrison and H.R. Quinn [BABAR Collaboration], The BaBar Physics Book, SLAC-R-0504, 
http://www.slac.stanford.edu/pubs/slacreports/slac-r-504.html. 
A. Ah, G. Kramer, and CD. Lu, Phys. Rev. D 58, 094009 (1998). 
J.A. Oiler and E. Oset, Phys. Rev. D 60, 074023 (1999). 
U.-G. Meifiner and J.A. Oiler, Nucl. Phys. A 679, 671 (2001). 
J.A. Oiler and E. Oset, Nucl. Phys. A 620, 438 (1997); ibid. 652, 407(E) (1997). 
O. Babelon, J.L. Basdevant, D. Caillerie, and G. Mennessier, Nucl. Phys. B 113, 445 (1976). 
J. Gasscr and H. Lcutwylcr, Annals Phys. 158, 142 (1984). 
J.F. Donoghuc, J. Gasscr, and H. Lcutwylcr, Nucl. Phys. B 343, 341 (1990). 
E.M. Aitala et al. [E791 CoUaboration], Phys. Rev. Lett. 86, 770 (2001). 
R. Casalbuoni et al., Phys. Rept. 281, 145 (1997); references therein. 
D. Becirevic and A. Le Yaouanc, J. High Energy Phys. 9903, 021 (1999). 

G. Burdman and J.F. Donoghue, Phys. Lett. B 280, 287 (1992); M.B. Wise, Phys. Rev. D 45, 2188 
(1992); T.-M. Yan et al, Phys. Rev. D 46, 1148 (1992); ibid. 55, 5851(E) (1992). 

A. F. Falk, Nucl. Phys. B 378, 79 (1992); A.F. Falk and M.E. Luke, Phys. Lett. B 292, 119 (1992); 
U. Kilian, J.G. Korner, and D. Pirjol, Phys. Lett. B 288, 360 (1992). 
L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

S. Ahmed et al. [CLEO Collaboration], Phys. Rev. Lett. 87, 251801 (2001). 
N.G. Deshpande, G. Eilam, X.-G. He, and J. Trampetic, Phys. Rev. D 52, 5354 (1995). 

B. Aubert et al. [BABAR Collaboration], arXiv:hcp-cx/0206004. 

H. -Y. Cheng and K.-C. Yang, arXiv:hep-ph/0205133. 



16 



